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SYSTEM AND METHOD FOR VISUALIZING CONDUCTIVITY 
AND CURRENT DENSITY DISTRIBUTION IN OBJECT 

5 Technical Field 

The present invention relates to system and method for visualization of 
conductivity (or resistivity) and/or current density distributions of a measuring object, 
such as a human body, or a substance. 

10 Background Art 

In general, X-ray, MRI, or ultrasound is used in visualization of an inside 
structure of the human body or the substance. However, the methods can not visualize 
electrical properties of the human body or the substance. 

As an effort for solving the problem, a method for visualizing current density 

15 distribution of an inside of a measuring object using the MRI is initially suggested by 
a research team of the Toronto University in 1989, and, thereafter, there have been 
active researches related thereto. This is the Current Density Imaging (CDI) technique. 
In the present CDI, a magnetic flux density B due to an injected current I is measured 
by using MRI technique, and a current density J is calculated by using the Ampere's 

20 law J = — VxB, for visualizing an inside current density distribution. 
Mo 

However, the CDI has a drawback in that rotation of the measuring object in 
an MRI scanner is required for obtaining the magnetic flux density B having three 
components Bx, By, and Bz. This is because the MRI scanner can only measure a z- 
directional component, a direction the same with a main magnetic field, of the flux 
25 density, i.e., Bz, when the measuring object is in the MRI scanner. That is, since the 
component the MRI scanner can measure at a time is only the Bz component, for 
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obtaining all the required three components of the magnetic flux density vector, the 
present CDI technique has a serious drawback of requiring the measuring object (the 
human body, or the substance) to be rotated in the MRI scanner. 

Meanwhile, as a known method for visualizing the electrical properties of the 
5 human body, or the substance, there has been Electrical Impedance Tomography (EIT) 
that has been under active research starting from late 1970s. The EIT provides an 
image of resistivity (or conductivity) distribution, an electrical property of the 
measuring object. 

In the EIT mostly taking the human body as the measuring object, many 
10 electrodes are attached on a surface of the human body, for visualizing the resistivity 
distribution of an inside of the human body. The visualization of the human body 
according to the resistivity is made possible since tissues of the human body, such as 
blood, bones, lung, heart, and the like have electrical properties different from one 
another. 

15 However, due to a fundamental drawback of the EIT, an image of the EIT is 

poor, to support only a low resolution. That is, EIT has a fundamental drawback in 
that a current-voltage data measured by EIT is extremely insensitive to variation of 
resistivity of an inside of the human body. Therefore, clinical application of EIT is not 
active, presently. 

20 

Disclosure of Invention 

An object of the present invention, designed to solve the related art problem, 
is to provide system and method for visualization of conductivity and current density 
distribution, in which electrical properties (conductivity and current density) of a 
25 measuring object can be visualized in more accurate high resolution without rotating 



WO 2005/051 163 PCT/KR2004/003026 

3 

the measuring object. 

Another object of the present invention is to provide system and method for 
visualization of conductivity and current density distribution, in which noises can be 
removed from magnetic flux density signals effectively and a high quality 
5 conductivity image can be produced. 

The object of the present invention can be achieved by providing a system for 
visualizing conductivity and current density distributions including a plurality of 
current injecting devices for injecting currents into a measuring object, a measuring 
unit for measuring a magnetic flux density due to the currents injected into a 
10 measuring object, an operating unit for selecting one pair of the current injecting 
devices in succession so as to inject currents of different directions into the measuring 
object, and calculating directional components of an anisotropic conductivity inside of 
the measuring object on the basis of the measured magnetic flux density, and 
displaying means for visualizing an inside of the measuring object by using the 
15 calculated directional components of the anisotropic conductivity. 

The selected pair of the current injecting devices apply the current to the 
measuring object, and one or more non-selected pairs of the current injecting devices 
measure a surface voltage of the measuring object. 

The measuring unit measures the magnetic flux density corresponding to one 
20 direction parallel to a main magnetic field thereof. 

The operating unit calculates an inside voltage and a surface voltage of the 
measuring object on the basis of an arbitrary conductivity, and calculates the 
directional components of the anisotropic conductivity on the basis of the calculated 
inside voltage and the measured magnetic flux density. 
25 The operating unit calculates new directional components of the anisotropic 
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conductivity on the basis of a calculated inside voltage depending on the calculated 
directional components, if an absolute value of a difference between the arbitrary 
conductivity and the calculated anisotropic conductivity is greater than a preset value. 

The operating unit determines that the calculated anisotropic conductivity is a 
5 true conductivity, if an absolute value of the difference between the arbitrary 
conductivity and the calculated anisotropic conductivity is smaller than a present 
value. 

The operating unit calculates a current density distribution inside of the 
measuring object on the basis of the calculated inside voltage and the measured 

10 magnetic flux density, and calculates a z-directional component on the basis of the 
calculated inside voltage and the calculated current density distribution. 

The operating unit divides the measured magnetic flux density into a first part 
not including noises and a second part including the noises, removes the noises from 
the second part, and combines the first part and the second part. 

15 In another aspect of the present invention, a method for visualizing 

conductivity and current density distributions includes (a) injecting currents of 
different directions into a measuring object through current injecting devices attached 
to a surface of the measuring object, (b) measuring a magnetic flux density due to 
the injected currents, (c) calculating directional components of an anisotropic 

20 conductivity inside of the measuring object on the basis of the measured magnetic flux 
density, and (d) visualizing an inside of the measuring object by using the calculated 
directional components of the anisotropic conductivity. 

The injecting currents of different directions into the measuring object 
comprises selecting a pair of the current injecting devices in succession, injecting the 

25 currents into the measuring object in succession through the selected pair of the 
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current injecting devices, and measuring a surface voltage of the measuring object 
through non-selected current injecting devices. 

The measuring the magnetic flux density comprises measuring the magnetic 
flux density corresponding to one direction parallel to a main magnetic field. 
5 The calculating the directional components of the anisotropic conductivity 

comprises calculating an inside voltage and a surface voltage of the measuring object 
on the basis of an arbitrary conductivity, and calculating the directional components of 
the anisotropic conductivity on the basis of the calculated inside voltage and the 
measured magnetic flux density. 

10 The calculating the directional components of the anisotropic conductivity 

comprises multiplying or dividing the calculated directional components of the 
anisotropic conductivity by a ratio of the calculated surface voltage and a measured 
surface voltage of the measuring object. 

The calculating the directional components of the anisotropic conductivity 

15 comprises calculating an absolute value of a difference of the arbitrary conductivity 
and the calculated anisotropic conductivity, and comparing the absolute value of the 
difference to a preset value. If the absolute value of the difference of the arbitrary 
conductivity and the measured anisotropic conductivity is greater than the preset value, 
new directional components of the anisotropic conductivity are calculated on the basis 

20 of a calculated inside voltage depending on the calculated directional components. If 
the absolute value of the difference of the arbitrary conductivity and the calculated 
anisotropic conductivity is smaller than a present value, it is determined that the 
calculated anisotropic conductivity is a true conductivity. 

The calculating directional components of the anisotropic conductivity 

25 comprises calculating five directional components of the anisotropic conductivity on 
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the basis of the calculated inside voltage and the measured magnetic flux density, 
calculating a current density distribution inside of the measuring object on the basis of 
the calculated inside voltage and the measured magnetic flux density, and calculating 
one directional component of the anisotropic conductivity on the basis of the 
5 calculated inside voltage and the calculated current density distribution. 

In another aspect of the present invention, a system for visualizing 
conductivity and current density distributions includes a plurality of current injecting 
devices for injecting currents into a measuring object, a measuring unit for measuring 
a magnetic flux density due to the currents injected into a measuring object, and an 

10 operating unit for dividing the measured magnetic flux density into a first part not 
including noises and a second part including the noises, and combining the first part 
and the second part after removing the noises from the second part. 

The operating unit removes the noises from the second part on the basis of 
an intensity variation of an MRI magnitude image, or removes the noises from the 

15 second part on the basis of a. spatial gradient of an image corresponding to the second 
part. 

The operating unit calculates a magnetic flux density of a portion in which the 
magnetic flux density is not measured, on the basis of the magnetic flux density 
measured at a border of the portion. 

20 In another aspect of the present invention, a method for visualizing 

conductivity and current density distributions comprises (a) injecting currents into a 
measuring object through current injecting devices attached to a surface of the 
measuring object, (b) measuring a magnetic flux density of due to the injected currents, 
(c) dividing the measured magnetic flux density into a first part not including noises 

25 and a second part including the noises, (d) removing the noises from the second part, 
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and (e) combining the first part and the second part. 

The dividing the measured magnetic flux density into the first part and the 
second part comprises calculating the first part by applying the measured magnetic 
flux density to a preset algorithm, and calculating a difference between the measured 
5 magnetic flux density and the calculated first part. 

The removing the noises from the second part comprises removing the noises 
from the second part on the basis of an intensity variation of an MRI magnitude image, 
or removing the noises from the second part on the basis of a spatial gradient of an 
image corresponding to the second part. 

10 

Brief Description of Drawings 

The accompanying drawings, which are included to provide a further 
understanding of the invention, illustrate embodiment(s) of the invention and together 
with the description serve to explain the principle of the invention. In the drawings; 
15 FIG 1 illustrates a diagram of a system for visualizing conductivity and 

current density distribution in accordance with a preferred embodiment of the present 
invention; 

FIG 2 illustrates a diagram showing a current injecting device in FIG 1 in 
accordance with a first preferred embodiment of the present invention; 
20 FIGS. 3 A and 3B illustrate diagrams each showing another embodiment of the 

current injecting device in FIG 1; 

FIGS. 4 illustrates a diagram showing each showing current flows in the 
measuring object and each flow corresponds to one of six different injection currents; 

FIG 5 illustrates a flow chart showing a method for visualization of 
25 conductivity and current density distribution in accordance with the present invention; 
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FIG 6A illustrates an ideal image of the magnetic flux density in which noises 
are not included; 

FIG 6B illustrates an image of the magnetic flux density in which noises are 
included; 

5 FIG 6C illustrates a second part of the magnetic flux density in which 

information representing delicate patterns and edges in the measuring object is 
included and noises are included; and 

FIG 6D illustrates a first part of the magnetic flux density in which baseline 
information representing sluggish patterns and edges in the measuring object is 
included. 

Best Mode for Carrying Out the Invention 

Reference will now be made in detail to the preferred embodiments of the 
present invention, examples of which are illustrated in the accompanying drawings. In 
describing the embodiments, same parts will be given the same names and reference 
symbols, and repetitive description of which will be omitted. FIG 1 illustrates a 
diagram of a system for visualizing conductivity and current density distribution in 
accordance with a preferred embodiment of the present invention. 

Referring to FIG 1, the system for visualizing conductivity and current 
density distributions includes current injecting devices 100 for injecting currents in 
directions different from one another into a measuring object, such as a human body 
or a substance, a magnetic flux density measuring unit (e.g., MRI scanner) 200 for 
measuring a magnetic flux density due to each of the currents, an operating unit 300 
for calculating a conductivity (or resistivity) distribution and a current density 
distribution inside of the measuring object from the currents and the magnetic flux 
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densities induced thereby, displaying means 400 for displaying the inside of the 
measuring object as an image of the conductivity and current density distributions, 
and a current source 500 for supplying the currents to the current injecting device 100, 
The current injecting device 100 includes an electrode 21 formed of a non- 
5 magnetic metal, such as copper, or a material having very high conductivity, and a 
wire for supplying the current to the electrode 21. Though direct attachment of the 
electrode 21 to a surface of the measuring object is possible, it is preferable that the 
electrode 21 is spaced apart from the surface of the measuring object for preventing 
distortion of the magnetic flux density at the surface of the measuring object. 

10 Therefore, as shown in FIG 2, one embodiment of the current injecting device 100 of 
the present invention further includes an insulating container 25 having the electrode 
21 attached thereto for putting the electrode 21 on a place spaced apart from the 
measuring object. The insulating container 25 has one surface having the electrode 21 
attached thereto, and the other surface to be brought into contact with the measuring 

15 object in an opened state. The insulating container 25 is stuffed with sponge having an 
electrolyte gel or electrolyte solution absorbed therein. The sponge having an 
electrolyte gel or electrolyte solution absorbed therein makes uniform flow of the 
current from the electrode 21 in the insulating container 25. Thus, by putting the 
electrode on a place spaced apart from the surface of the measuring object by using 

20 the insulating container 25, and supplying the current to the measuring object 
uniformly by using the sponge having the electrolyte gel or electrolyte solution 
absorbed therein, distortion of an MRI image at the surface of the measuring object in 
contact with the electrode 21 can be reduced. 

Referring to FIGS. 3A and 3B, another embodiment of the current injecting 

25 device 100 of the present invention includes an electrode 31, a wire 33 connected to 
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the electrode 31, an semispherical insulating container 35 having the electrode 31 
attached thereto. The insulating container 35 has one side having the electrode 31 
attached thereto, and the other surface to be brought into contact with the measuring 
object 'S' in an opened state. The insulating container 35 is stuffed with sponge 37 
5 having the electrolyte gel or electrolyte solution absorbed therein. The sponge 37 may 
be replaced by only the electrolyte gel. 

Referring to FIG 4, the current injecting devices 100 are attached to a 
circumference of the measuring object C S\ and the current is made to flow only to one 
pair of the current injecting devices 100 at a time, such that the current flows through 

10 the inside of the measuring object 'S' with the one pair of the current injecting devices 
100. If it is assumed that a number of the current injecting devices 100 is 4 E\ a 
number 'N' of cases of the pairs of the current injecting devices 100 for injecting the 
current into the measuring object *S' is E(E-l)/2. This is equal to a number of current 
paths in the measuring object *S\ A current between the pair of current injecting 

15 devices 100 is represented with V (j = 1, 2, — , N). The current V induces a current 
density J* (J*x, J*y, fz) in the measuring object *S\ 

The MRI scanner 200 measures a magnetic flux density B* in the measuring 
object 'S' due to the current I j , particularly, a z-directional magnetic flux density B?z 
parallel to a main magnetic field of the MRI scanner 200. That is, the MRI scanner 

20 200 measures magnetic flux densities B l z, B 2 z, B 3 z, — , B N z due to currents I 1 , 1 2 , 1 3 , 
— , I N injected into the measuring object 6 S\ respectively. 

The operating unit 300 selects one of the pairs of the current injecting devices 
100 in succession, and controls the one of the pair of the current injecting devices to 
inject the current into the measuring object 6 S\ The operating unit 300 also applies the 

25 magnetic flux density B*z measured by the MRI scanner 200 to a preset algorithm, and 
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calculates a conductivity a and a current density J inside of the measuring object 'S\ 
Since the conductivity a is an inversion of the resistivity p, the conductivity a can be 
represented with 1/p. Then, the displaying means 400 visualizes the inside of the 
measuring object 'S' according to the conductivity a and the current density JP. 
5 Optionally, the operating unit 300 may remove the noises from the measured 

magnetic flux density B*z. The operating unit 300 divides the measured magnetic flux 
density B*z into a first magnetic flux density part not including the noises and a second 
magnetic flux density part including the noises, and synthesizes the first magnetic flux 
density part and the second magnetic flux density part after removing the noises from 
10 the second magnetic flux density part. 

A method for visualizing conductivity and current density distributions in 
accordance with a preferred embodiment of the present invention, and an algorithm 
applied thereto will be described, in detail. 



15 FIRST EMBODIMENT 

The current injecting devices 100 are attached to a circumference of a 
measuring object 'S', i.e., a human body, and injects the currents in six directions into 
the measuring object 'S' in order to induce a three-dimensional anisotropic 
conductivity. The currents to be injected in six directions are generated by at least the 

20 four current injecting devices. A matrix of the anisotropic conductivity A is as follows. 

a U a 12 a 13 
a 21 a 22 CT 23 
CT 31 a 32 a 33. 

The anisotropic conductivity 'A' has nine directional components. The 
anisotropic conductivity 'A' practically includes the six directional components 
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because the component <j i2 is identical with the component ^- 2| , the component 

is identical with the component and the component ^is identical with 

the component ^ n . Therefore, the current injecting devices 100 inject the currents in 

at least the six directions into the measuring object 'S' in order to induce the 
5 anisotropic conductivity including the six components. A method for visualizing 
conductivity and current density distributions in accordance with an example in which 
the currents are injected into the measuring object *S' in seven directions, will be 
described. 

The measuring object *S' is placed in an MRI scanner 200, and the pairs of 
10 the current injecting devices 100 are selected one by one in succession, and the current 
I J (j = 1, 2, — , 7) is supplied to an inside of the measuring object 6 S* through the pair 
of the current injecting devices 100 selected in succession (Sll). For an example, as 
shown in FIG 4, in a case in which four current injecting devices 100a, 100b, 100c, 
and lOOd are attached to a circumference of the measuring object 6 S\ the currents are 
15 supplied to the inside of the measuring object 'S' in six directions in succession.. In 
these instances, the MRI scanner 200 measures the magnetic flux densities B x z, 
B 2 z,... B 7 z due to the currents I 1 , I 2 ,... I 7 injected into the measuring object 'S' 
respectively. The MRI scanner 200 measures 4 z' direction components of the magnetic 
flux densities, only. Also, the non-selected current injecting devices 100 measure a 
20 surface voltage V j |as of the measuring object 'S' (S12). 

Then, the operating unit 300 applies the magnetic flux densities B*z, 
B 2 z,...B 7 z measured by the MRI scanner 200 and the surface voltage V j |as measured 
by the non-selected current injecting devices to a preset algorithm, and calculates the 
anisotropic conductivity 'A* and the current density J inside of the measuring object 
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Equations to be used for the algorithm in the first embodiment will be 
described, where, 

S : a measuring object, 

9S : a surface of the measuring object, 

g 1 : a component of a current density of an injected current at a surface of the 
measuring object, 

n : a unit normal vector at the surface of the measuring object, 

V J : a voltage of the injected current in the measuring object. 

The magnetic flux densities B l z, B 2 z,. . .B 7 z can be expressed as the following 

equation (1) based on the Ampere's law - — V7 2 £ = Vxj and the Ohm's law 

Mo 

J = -AW . The symbol is a magnetic permeability in a free space. 



1 

?0 



V 2 B] 









d x v* d y v> 9,v* dy-d 2 y v* d ,d ^ - d y a y 
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-(1) 



If the measured magnetic flux density B J z and the voltage V J are substituted 
into the equation (1), the components of the anisotropic conductivity 'A' may be 
obtained. However, it is impossible to obtain the actual components of the anisotropic 
conductivity 6 A' because the voltage V j inside the measuring object *S' cannot be 
measured. In order to obtain the voltage V j , an iterative algorithm is used. If the 
20 voltage values are obtained by using the iterative algorithm, the components ff , 

<Jn> CTu' CTn 9 °* * e anisotropic conductivity 4 A' except the component 
may be calculated by substituting the voltage values into the equation (1). 
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The component Q- may be calculated by substituting the calculated 
components q- 9 <j i2 , q-^, ^ , the voltage V j , and the current densities J*x, 

J*y into the following equation (2). 



(2) 





d 2 y 










( * £a) 













Wherein, the current densities Jx, J*y are may obtained from the following 
equations. 



o n a n a 13 

a 21 a 22 a 23 
a 3l a 32 °33. 



K*\ 



W 

10 The value of the voltage V J should be obtained in order to calculate the 

components of the anisotropic conductivity by using the equations (1), (2). The value 
of the voltage may be obtained by using an iterative algorithm and the following 
equation (3). 

v • (A v^) = o (/=U,...7) 

15 -AV^n^ (3) 

The equation (3) is to calculate the voltage V j inside the measuring object 'S\ 
The component g* of the current density is zero (0) at the surface of the 
measuring object 'S* not in contact with the selected pair of the current injecting 
devices 100. If the anisotropic conductivity 6 A\ the injected current P, and positions 
20 of the current injecting devices 100 are known, the voltage V j can be calculated by a 
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calculating method, such as FEM (Finite Element Method). 

The anisotropic conductivity 'A' should be known in order to calculate the 
voltage V j by using the equation (3), The voltage V j should be known in order to 
calculate the anisotropic conductivity 'A' by using the equations (1), (2). A process for 
5 calculating the conductivity 'A' and the voltage V j at the operating unit 300 by using 
the foregoing equations (1), (2), (3) will be described below in detail. 

The initial conductivity A m (m=0) is set, and the initial conductivity A m (m=0) 
and the current density component g* at the surface of the measuring object 6 S' are 
substituted into the following equation (4), to calculate voltages V^i and V j m +i|as at 
10 the inside and the surface of the measuring object *S' (S13). The initial conductivity 
Am (m=0) may be a constant value if it is assumed that the initial conductivity Am is 
isotropic, or may has the directional components acknowledged according to the 
tissues of the body. The inner voltage VVi and the surface voltage Villas are 
calculated based on the initial conductivity Am. The current density component g* can 
15 be known from the current F (j = 1 , 2, — , 7). 

-A^vUts-n-g* (4) 

After the inner voltage V j m +i corresponding to the initial conductivity is 
obtained, the inner voltage V j m +i and the measured magnetic flux density & z are 
20 substituted into the equation (1) and the equation (2) in order to calculate the 
components of the anisotropic conductivity A m +i (SI 4). The calculated anisotropic 
conductivity A m +i is closer to a true conductivity than the conductivity A m , because 
the calculated conductivity Am+i has taken an true magnetic flux density B j 2 inside of 
the measuring object *S' into account. 
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Moreover, in the present invention, for obtaining the conductivity close to the 
true value, an iterative form of algorithm is used. That is, alike the following equation 
||cr OT+1 -o-J^s , a difference between a conductivity A ra in a prior step and a 

conductivity Am+i in the present step is calculated, and it is determined whether an 
5 absolute value of the difference is smaller than a preset value c e' or not (SI 5). The 
equations (1), (2), and (4) are calculated repeatedly by increasing the *m' until the 
calculated difference becomes smaller than the preset value *e' (SI 6). If the absolute 
value of the difference of the conductivity A m in the prior step and the conductivity 
Am+i in the present step is smaller than the preset value *e\ the conductivity Am+i is 

10 taken as the true conductivity. Then the measuring object 6 S' is visualized by using the 
calculated conductivity A m +i (SI 8). 

However, even if the calculated conductivity Am+i has taken the true magnetic 
flux density B* z into account, the calculated conductivity Am+i may differ from the true 
conductivity. That is, though the calculated conductivity A m +i has a value varied with 

15 the magnetic flux density B^, the calculated conductivity A m +i may be greater or 
smaller than the true conductivity by a multiplication of a constant. For eliminating 
the difference between the calculated conductivity Am+i and the true conductivity, in 
the present invention, the measured surface voltage V j |as and the calculated surface 
voltage V m+ i j |as may be used for scaling the conductivity A m +i. In other words, 

20 according to a ratio of the measured surface voltage V J |as to the calculated surface 
voltage Vm+i^as, the conductivity A m +i is multiplied or divided by a constant value. 
Thus, the conductivity Am+i is calculated correctly. 

In addition, the conductivity A^i is substituted into the equation (4) to 
calculated a voltage depending on the conductivity Am+i, i.e., the voltage V j inside of 
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the measuring object *S' is calculated. Then, by using the conductivity Am+i and the 
voltage V j , the current density J* is calculated, and by using the current density J*, the 
inside of the measuring object 'S 5 is visualized. 



5 SECOND EMBODIMENT 

The current is injected into the measuring object 6 S' in two directions or more 
in succession in order to induce a conductivity a, and is generated by at least 3 current 
injecting devices 100. 

The measuring object 4 S' is placed in the MRI scanner 200, and the pairs of 
10 the current injecting devices 100 are selected one by one in succession, and the current 

V (j = 1, 2, — , N) is supplied to an inside of the measuring object 'S' through the pair 
of the current injecting devices 100 selected in succession. The MRI scanner 200 
measures the magnetic flux densities B*z, B 2 z,... B N z due to the currents I 1 , I 2 ,...I N 
injected into the measuring object 'S' respectively, the non-selected current injecting 

15 devices 100 measure a surface voltage V j |as of the measuring object 4 S' when the 
selected pair of the current injecting devices 100 injects the current into the measuring 
object 'S\ 

Then, the operating unit 300 applies the magnetic flux densities B*z, 
B 2 z,...B N z measured by the MRI scanner 200 and the surface voltage V j |$s measured 
20 by the non-selected current injecting devices to a preset algorithm, and calculates the 
conductivity a and the current density J 5 inside of the measuring object *S*. 

The conductivity a is calculated on the basis of the Maxwell equation 

V x b? = n 0 j J including a magnetic flux density and a current density. The Maxwell 
equation is multiplied by a vector function * , and is integrated three times, then the 
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following equation (5) depending on the Ohm's law f—Q^ V' mav t, e obtained. 

Js °Js (5) 

Wherein, ^foj^) . Since the only component B j z can be measured 
without rotating the measuring object *S\ the components B'x and B*y of the magnetic 
5 flux density should be removed from the equation (5). Thus, the function group, such 
as the following equation (6), is determined. 



T:= 



•y 

4^"(r)=0 ifnxT=0 
dz 3 



V G H\S): foraU^G d S, ^(Vo^Wif n • "e!=0 

a x ay 



(6) 

Wherein, ^a 3 * 0 * 0 ' 1 ) . ^(S) i s a square integrable function. The 
divergence theorem induces the following equation (7) by substituting the function 

10 dy ' dx 9 into the equation (5). 



<7) 



Wherein, the symbol V i s dx' dy and the symbol L is l-i OJ . 

If the measured magnetic flux density B*z and the voltage V J are substituted 
into the equation (7), the conductivity a may be obtained. However, it is impossible to 
15 obtain the actual conductivity because the voltage V j inside the measuring object 'S' 
cannot be measured. In order to obtain the voltage V j , an iterative algorithm based on 
the following equation (8) is used. 



-(8) 



Wherein, the conductivity a m +i is the newly calculated conductivity and the 
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voltage V j m is the calculated voltage depending on the previous conductivity a m . The 
voltage V j ra may be calculated by the following equation (9) and a calculating method, 
such as FEM (Finite Element Method). 

V, K V W inside S 

5 -° W VK;. M=g* 9 tmdS 

J d s (9) 

The component g* of the current density is zero (0) at the surface of the 
measuring object 'S' not in contact with the selected pair of the current injecting 
devices 100. 

10 A process for calculating the conductivity a^i and the voltage V m at the 

operating unit 300 by using the foregoing equations (8) and (9) will be described 
below in detail. 

The initial conductivity a m (m=0) is set, and the initial conductivity a m and 
the current density component g* at the surface of the measuring object 'S' are 

15 substituted into the following equation (9), to calculate voltages V j m and V j m |as at the 
inside and the surface of the measuring object 4 S\ The initial conductivity a m may be 
a constant value if it is assumed that the initial conductivity o m is isotropic, or may has 
the directional components acknowledged according to the tissues of the body. The 
inner voltage V j m and the surface voltage V^las are calculated based on the initial 

20 conductivity a m . The current density component g* can be known from the current V (j 
= 1,2, — ,N). 

After the inner voltage V j m corresponding to the initial conductivity a m is 
obtained, the inner voltage Vm and the measured magnetic flux density B j z are 
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substituted into the equation (8) in order to calculate the conductivity a m +i. The 
calculated am+i is closer to a true conductivity than the conductivity a m , because the 
calculated conductivity Cm+i has taken an true magnetic flux density B* z inside of the 
measuring object 'S' into account 
5 Moreover, in the present invention, for obtaining the conductivity close to the 

true value, an iterative form of algorithm is used. That is, alike the following equation 

f f f \°m+x(h*Jry( 2 di 

J J J s , a difference between the conductivity a m in a prior step and 

the conductivity Gm+i calculated in the present step is calculated, and it is determined 
whether an absolute value of the difference is smaller than a preset value 4 e' or not. 

10 The equations (8) and (9) are calculated repeatedly by increasing the 'm' until the 
calculated difference becomes smaller than the preset value *e\ If the absolute value 
of the difference of the conductivity a m in the prior step and the conductivity am+i in 
the present step is smaller than the preset value *e\ the conductivity a m +i is taken as 
the true conductivity. Then the measuring object 'S' is visualized by using the 

1 5 calculated conductivity a^i . 

However, even if the calculated conductivity a m +i has taken the true magnetic 
flux density B* z into account, the calculated conductivity <Tm+i may differ from the true 
conductivity. For example, though the calculated conductivity a m +i has a value varied 
with the magnetic flux density B^z, the calculated conductivity Cm+i may be greater or 

20 smaller than the true conductivity. For eliminating the difference between the 
calculated conductivity Om+i and the true conductivity, in the present invention, the 
measured surface voltage V j |as and the calculated surface voltage VnnVlas ma V ^ e usec ^ 
for scaling the conductivity a^i. In other words, the conductivity am+i is multiplied 
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or divided by a ratio of the measured surface voltage V j |as to the calculated surface 
voltage V m +i J |as- Thus, the conductivity CTm+i is calculated correctly. 

In addition, the conductivity Om+i is substituted into the equation (9) to 
calculated a voltage depending on the conductivity a,n+i, i.e., the voltage V j inside of 
5 the measuring object 'S' is calculated. Then, by using the conductivity Om+i and the 
voltage V j , the current density J is calculated, and by using the current density J 1 , the 
inside of the measuring object 'S' is visualized. 

A process for calculating the conductivity based on first partial derivatives of 
the measured magnetic flux density B j z using another equations will be described 
10 below in detail. 

The vector * °) can be expressed as the following equation (10). 
(-4,^,0)= Vw> + VX0> m (1Q) 

Wherein, the symbol ^ is a scalar function and the symbol ^ is a 

WCr) = f f — l^r- -y-[aVJ^(P)]<£4 
vector function J J ss 4n\r-r] dz . The symbol is 

15 referred to a square measure. 

The following equation (11) representing the relation between the flux density 

& z and the scalar function ^ may be obtained by using the Maxwell equation 
= y i nc i u ding the magnetic flux density ^ =ss (^^'^) and the 
current density "^"(^'^^O . 

20 ^ 



WO 2005/051163 PCT/KR2004/003026 

22 

V 2 A / = 0 in 5 



r • Vtf = («„«„ 0) • J T • on 5 5 

VV = 0 in 5 



t v*V = U 



r on ds 



-(11) 



Wherein, the symbol T is a tangent vector ^ ( _, V w z>0) corresponding 

to 3S . 

The following equation (12) may be obtained by substituting the equation (11) 
into the equation (10). 

j{=M-+ L m , j^^-hv ] 

Wherein, the symbol ^ represents ^° , the function L lV\- 



represents dy & x d z , and the function M**] : represents 

10 dx dy dz , xhe functions and MKl: are determined on the 

basis of the measured data V j , g" , and & z . 

The following equation (13) may be obtained by using the equation (12) and 
the Ohm's law. 

dV> I dB 1 ^rrrJ^ 9 V I d if 



a — 



(13) 

15 If the measured magnetic flux density &z and the voltage V j are substituted 

into the equation (13), the conductivity a may be obtained. However, it is impossible 
to obtain the actual conductivity because the voltage V j inside the measuring object 'S' 
cannot be measured. In order to obtain the voltage V j , an iterative algorithm based on 
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the following equation (14) is used. 




... (14) 



Wherein, the conductivity a m+ i is the newly calculated conductivity and the 



voltage V J m is the calculated voltage depending on the previous conductivity a a The 



5 voltage V j m may be calculated by using a calculating method depending on the 
following equation (1 5), such as FEM (Finite Element Method). 

V«(a OT VKg=0 inside S 
■oX 1 ^ ondS 



measuring object 'S' not in contact with the selected pair of the current injecting 
devices 100. 

A process for calculating the conductivity Cm+i and the voltage V j m at the 
operating unit 300 by using the foregoing equations (14) and (15) will be described 
15 below in detail. 

The initial conductivity a m (m=0) is set, and the initial conductivity a m and 
the current density component g i at the surface of the measuring object 4 S* are 
substituted into the following equation (15), to calculate voltages V j ra and V j m |as at the 
inside and the surface of the measuring object *S\ The initial conductivity a m may be 
20 a constant value if it is assumed that the initial conductivity a m is isotropic, or may has 
the directional components acknowledged according to the tissues of the body. The 




(15) 



10 



The component g* of the current density is zero (0) at the surface of the 
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inner voltage V j m and the surface voltage V j m |as are calculated based on the initial 
conductivity a m . 

After the inner voltage corresponding to the initial conductivity a m is 



5 substituted into the equation (14) in order to calculate the conductivity cwi. The 
calculated a m +i is closer to a true conductivity than the conductivity a m , because the 
calculated conductivity a^H has taken an true magnetic flux density B* z inside of the 
measuring object C S' into account. 



10 true value, an iterative form of algorithm is used. That is, alike the following equation 



the conductivity am+i calculated in the present step is calculated, and it is determined 
whether an absolute value of the difference is smaller than a preset value 'e' or not. 
The equations (14) and (15) are calculated repeatedly by increasing the 'm' until the 

15 calculated difference becomes smaller than the preset value c e\ If the absolute value 
of the difference of the conductivity a m in the prior step and the conductivity a m+ i in 
the present step is smaller than the preset value 6 s', the conductivity a m +i is taken as 
the true conductivity. Then the measuring object 4 S' is visualized by using the 
calculated conductivity a ra+ i. 

20 However, even if the calculated conductivity cwi has taken the true magnetic 

flux density & z into account, the calculated conductivity a^i may have scaling 
difference from the true conductivity. For example, though the calculated conductivity 
Ojn+x has a value varied with the magnetic flux density B J Z , the calculated conductivity 



obtained, the inner voltage V j m and the measured magnetic flux density B j z are 



Moreover, in the present invention, for obtaining the conductivity close to the 




, a difference between the conductivity a m in a prior step and 
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a m +i may be greater or smaller than the true conductivity. For eliminating the scaling 
difference between the calculated conductivity Qm+i and the true conductivity, in the 
present invention, the measured surface voltage V j |as and the calculated surface 
voltage V m+ i j |as may be used for scaling the conductivity am+i. In other words, the 
5 conductivity a m +i is multiplied or divided by a ratio of the measured surface voltage 
V J |as to the calculated surface voltage Vm+^las. Thus, the conductivity a m+ i is 
calculated correctly 

In addition, the conductivity o^i is substituted into the equation (15) to 
calculated a voltage depending on the conductivity Gm+i, i.e., the voltage V inside of 
10 the measuring object 'S' is calculated. Then, by using the conductivity a^i and the 
voltage V j , the current density J* is calculated, and by using the current density J*, the 
inside of the measuring object 'S' is visualized. 

THIRD EMBODIMENT 

15 The measuring object 'S' is placed in the MRI scanner 200, and the pairs of 

the current injecting devices 100 are selected one by one in succession, and the current 
I J (j = 1, 2, — , N) is supplied to an inside of the measuring object 'S' through the pair 
of the current injecting devices 100 selected in succession. The MRI scanner 200 
measures the magnetic flux densities B*z, B 2 z,... B N z due to the currents I 1 , I 2 ,...I N 

20 injected into the measuring object 'S' respectively, the non-selected current injecting 
devices 100 measure a surface voltage V J |as of the measuring object 'S' when the 
selected pair of the current injecting devices 100 injects the current into the measuring 
object »S\ 

FIG 6A illustrates an ideal image of the magnetic flux density not including 
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the noises, FIG 6B illustrates an image of the magnetic flux density including the 
noises. Referring to FIG 6B, the measured magnetic flux density includes the 
noises due to the mechanical errors of the MRI scanner 200 and the environmental 
factors. In order to remove the noises from the measured magnetic flux density B' z , the 
5 operating unit 300 divides the measured magnetic flux density B j z into the first part 
B h z and the second part B f z . 

The operating unit 300 calculates the first part B h z of the measured magnetic 
flux density based on the following equation (16) satisfying the boundary condition in 
order to separate the first part B h z and the second part B f z . 

10 V2 *?=° (16) 

The operating unit 300 calculates the difference between the measured 
magnetic flux density B j z and the first part B h z in order to obtain the second part B^ 

FIG 6C illustrates the second part B f z of the magnetic flux density including 
information representing delicate patterns and edges in the measuring object and the 
15 noises, and FIG 6D illustrates the first part B h z of the magnetic flux density including 
baseline information representing sluggish the patterns and the edges in the measuring 
object. 

The noises should be removed from the second part B z of the magnetic flux 
density. In the present invention, for removing the noises from the second part B f z , two 
20 diffusion coefficients are suggested. A first diffusion coefficient is based on 
information related to the preliminary MRI magnitude image provided by the MRI 
scanner 200. The information related to the preliminary MRI magnitude image 

depends on the intensity variation ^ M of the preliminary MRI magnitude image, 
and includes the information representing the patterns and the edges in the measuring 
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object *S\ A second diffusion coefficient is based on the space gradient ^ u of the 
image corresponding to the second part B f z . The following equation (17) shows the 
diffusion model including the first and the second diffusion coefficients 

^-V«|VA4)1(|Vi©Vii] 

5 Wherein, an initial condition includes ^ x >y>ty~ Bf z( x >y) 9 and the function 

1 

*(V") represents Tjv^TT . 

The operating unit 300 may remove the noises from the second part B f z by 
using at least one of the two diffusion coefficients. The first diffusion coefficient is 
used to prevent from blurring B f z along the boundary of the tissues in the measuring 

10 object 'S' using the intensity variation of the preliminary MRI magnitude 

image. The second diffusion coefficient is used to preserve edge information of the 

image B f z using the space gradient ^ u of the image corresponding to the second 
part B f z . 

After the noises are removed from the second part B f z , the operating unit 300 
15 combines the first part B h z and the second part B f z . Thus the denoised magnetic flux 
density may be obtained. 

The magnetic flux density can not be measured in a portion of the measuring 
object *S' filled with the air. Therefore, the operating unit 300 detects the border 
between the portion with the air and the portion without the air from the MR 
20 magnitude image, and calculates the magnetic flux density B b z in the portions filled 
with the air on the basis of the magnetic flux density B^ measured at the border of the 
portions. That is, the operating unit 300 may calculate the magnetic flux density in the 
portion with the air by substituting the solution of the Laplace's equation 
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* in the portion with the boundary condition having the magnetic flux density 
measured at the border. Then the operating unit 300 applies the denoised magnetic 
flux density B j 2 and the calculated magnetic flux density B b 2 to the algorithm 
suggested in the first embodiment or the second embodiment in order to calculate the 
5 conductivity a and the current density J 1 . 



Industrial Applicability 

As has been described, the present invention permits to obtain a conductivity 
distribution image and a current density distribution image of an inside of a measuring 

10 object without rotating the measuring object, such as a human body, or a substance, 
within an MRI scanner. 

Moreover, the present invention permits to produce a high quality image by 
removing the noises from the measured magnetic flux density successfully and 
calculating the magnetic flux density in the tissue filled with the air or water on the 

15 basis of the magnetic flux density measured at the border of the tissue. 

Furthermore, the conductivity distribution image and the current density 
distribution image can be applicable to the present medical equipment. According to 
this, a functional imaging can be provided for organs each of which biological tissue 
conductivity varies with physiologic action, to secure a new medical examination 

20 technology. For an example, an image required for measuring a volume of a lung 
varies with respiration (a measurement by using variation of a resistivity of the lung in 
inhale and exhale of air), a measurement of an cardiac output of a heart, a 
measurement of a stomach, a measurement of cystic residual urine, a measurement of 
a healed state of a bone fracture, or a measurement of a brain function. 
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Since quantitative electric properties of a human tissue can be known, 
optimization of diagnosis equipments, and treatment equipments using electro- 
magnetic energy, such as various electric/magnetic stimulators, defibrillators, cardiac 
pacemakers, and the like, is made possible. 
5 It will be apparent to those skilled in the art that various modifications and 

variations can be made in the present invention without departing from the spirit or 
scope of the invention. Thus, it is intended that the present invention cover the 
modifications and variations of this invention provided they come within the scope of 
the appended claims and their equivalents. 



